Nonadiabatic extension of the Heisenberg model 
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The localized states within the Heisenberg model of magnetism should be represented by best 
localized Wannier functions forming a unitary transformation of the Bloch functions of the narrow- 
est partly filled energy bands in the metals. However, as a consequence of degeneracies between 
the energy bands near the Fermi level, in any metal these Wannier functions cannot be chosen 
symmetry-adapted to the complete paramagnetic group M p . Therefore, it is proposed to use Wan- 
nier functions with the reduced symmetry of a magnetic subgroup M of M [case (a)] or spin 
dependent Wannier functions [case (b)]. The original Heisenberg model is reinterpreted in order 
to understand the pronounced symmetry of these Wannier functions. While the original model 
assumes that there is exactly one electron at each atom, the extended model postulates that in 
narrow bands there are as many as possible atoms occupied by exactly one electron. However, this 
state with the highest possible atomiclike character cannot be described within the adiabatic (or 
Born-Oppenheimer) approximation because it requires a more realistic description of the electronic 
motion. Within the (true) nonadiabatic system the electrons move on localized orbitals that are 
still symmetric on the average of time, but not at any moment. These nonadiabatic states have the 
same symmetry as the adiabatic states and determine the commutation properties of the nonadi- 
abatic Hamiltonian H n . The nonadiabatic Heisenberg model is a purely group-theoretical model 
which interprets the commutation properties of H n that are explicitly given in this paper for the 
two important cases (a) and (b). There is evidence that the occurrence of these two types of Wan- 
nier functions in the band structure of a metal is connected with the occurrence of magnetism and 
superconductivity, respectively. 



PACS numbers: Heisenberg model, 75.10.-b, 74.20.-z 

I. INTRODUCTION 

In any application of the Heisenberg model of 
magnetism, 1 the localized states of the electrons should 
be represented by Wannier functions 
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with the following properties: 

(i) the u>i(r — R~ pi) are centered on the atomic posi- 
tions R + pi; 

(ii) the Wi(f — R— pi) are gained by a unitary trans- 
formation from the Bloch functions (ft (r) of the energy 
bands of interest; 

(hi) the Wi(f — R — pi) are symmetry- adapted to the 
space group G of the considered metal; and 

(iv) the Wi(r — R — pi) are as well localized as possible. 

The first sum in Eq. (1.1) runs over the N vectors k of 
the first Brillouin zone (BZ), the second sum runs over 
the fj, bands of interest (with the band indices q = 1 to p), 
and R and pi denote the vectors of the Bravais lattice and 
the centers of symmetry of the Wannier functions within 
the unit cell, respectively. 

The transformation is unitary [point (ii)] if the coef- 
ficients giq(k) in Eq. (1.1) are the elements of a unitary 
matrix g(k), 



The Wannier functions are symmetry-adapted to G 
[point (iii)] if they satisfy the equation 

Wi(a _1 (f- R- pi)) = y^^Dji(a)wj(r- R-pi) (1.3) 

for the elements a of the point group Go of G, where the 
matrices [Dji(a)] form a (reducible or irreducible) single- 
valued representation D of Go, see Eq. (1.8) of Ref. 2. 
[Note that in Eq. (1.3) on the right hand side there is pi 
and not pj.] 

The Wannier functions are as well localizable as possi- 
ble [point (iv)] if the Bloch-like functions 



9=1 



(1.4) 



- 1 {k)=g\k). 



(1.2) 



vary (for fixed r) smoothly through the whole k space. 2 
It is one of the most important results of the group the- 
ory of Wannier functions that Wannier functions comply- 
ing with all the four conditions given above exist only in 
isolated sets of p energy bands which satisfy the group- 
theoretical compatibility relations throughout the Bril- 
louin zone. 2 ' 3,4 ' 5,6 [This condition is necessary, but not 
sufficient: in addition, there must exist unitary matrices 
S(K) defined in Eq. (4.16) of Ref. 2 which satisfy the 
equations (4.17) and (4.28) of Ref. 2. These matrices 
S(K) determine the positions pi of the Wannier func- 
tions.] In an "isolated" set of energy bands, each band 
may be connected by degeneracies to the other bands of 



this set, but must not be connected to bands not belong- 
ing to this set. 

The mentioned "energy bands of interest" are the 
partly filled energy bands in the band structures of the 
considered metal. Often, it is only one roughly half-filled 
band which interests. However, the energy bands in the 
(paramagnetic) band structures of the metals are degen- 
erate at several points and lines of symmetry of the Bril- 
louin zone. Because of these degeneracies, is not possible 
to separate narrow isolated sets of energy bands which 
satisfy the compatibility relations throughout the Bril- 
louin zone. For this reason, Wannicr functions with all 
the properties demanded above do not exist in the metals. 

Therefore, the localized states are often represented by 
"approximated" Wannier functions which no longer form 
an exactly unitary transformation of the Bloch functions. 
These Wannier functions are constructed from slightly 
modified energy bands in which some of the Bloch func- 
tions at points, lines, and planes of symmetry arc re- 
placed by Bloch functions with a symmetry appropri- 
ate for the construction of Wannier functions. Hence, 
these approximated Wannier functions have lost all the 
information connected with the symmetry of the removed 
Bloch functions and carry the wrong information of the 
new Bloch functions. 

The nonadiabatic Hciscnbcrg model (NHM) as pro- 
posed in this and previous papers 7 ' 8 ' 9 ' 10 extends the orig- 
inal Hciscnbcrg model on the basis of Wannier func- 
tions which form an exactly unitary transformation of 
the Bloch functions of the bands of interest. Within this 
model, it is not allowed to replace any Bloch function in 
the calculated band structure by functions with a new 
symmetry. Hence, the NHM takes into account the com- 
plete information connected with the symmetry of the 
Bloch functions in the band structure of the considered 
metal. 

Clearly, the Wannier functions used within the NHM 
cannot not comply with all the properties (i) - (iv) given 
above. The development of the nonadiabatic model was 
suggested by two observations: 

1. An exactly unitary transformation of the Bloch 
functions of the partly filled bands into best localized 
Wannier functions becomes possible in nearly all the met- 
als when the Wannier functions are allowed to have the 
reduced symmetry of a magnetic subgroup M of the para- 
magnetic group [see appendix A, case (a)] or when they 
are allowed to be spin dependent [see appendix A, case 
(b)]- 

2. There is a deep connection between the symmetry of 
these exact Wannier functions in a metal and the physical 
properties of the electrons at the Fermi level. 7 ' 8 ' 9 ' 10 ' 11 

The original Hciscnbcrg model of magnetism is defined 
by the assumption that there is exactly one electron on 
each atom of a metal. The NHM replaces this assumption 
by introducing three postulates which will be given in 
the following Sec. II. These postulates combine in a new 
way the Heisenberg model with the band model. The 
fundamental second postulate given in Eq. (2.18) states 



that the probability to find exactly one electron on an 
atom is as large as possible in narrow energy bands. 

The second postulate of the NHM cannot be satisfied 
within the adiabatic (or Born-Oppcnhcimcr) approxima- 
tion. In the framework of this approximation the elec- 
trons move in rigid orbitals in the average potential of the 
other electrons. The second postulate, however, requires 
a more realistic description of the electronic motion. In 
the true (nonadiabatic) system a localized electron moves 
in a potential depending on which of the adjacent local- 
ized states is occupied and on the present motion of the 
electrons in these states. These modified orbitals 

\T,m,v) 

are described by introducing a new quantum number v 
which labels different states of motion of the center of 
mass of the localized states. 12 [T denotes the positions of 
the atoms, see Eq. (2.2).] 

Nonadiabatic localized functions 

(r,t,q\f,m,v) 

(as introduced in the next section) which represent the 
nonadiabatic localized states arc highly complicated. 
Hence, it will be practically impossible to give these func- 
tions explicitly. Fortunately, one important feature of 
these functions is known exactly: they have the same 
symmetry as the exact Wannier functions of the nar- 
rowest, roughly half- filled energy bands of the metal un- 
der consideration. Thus, any application of the NHM 
starts with a group-theoretical examination of the sym- 
metry of the best localized (spin dependent) Wannicr 
functions which is clearly determined by the symmetry 
of the Bloch functions in the band structure of the given 
metal. 2 ' 3 ' 4 ' 5 ' 6 ' 13 The symmetry of these Wannier func- 
tions is explicitly given for two important cases (a) and 
(b) in appendix A. 

The NHM is a purely group-theoretical model. An 
explicit knowledge of the nonadiabatic localized functions 
(going beyond of their symmetry) does not provide new 
physical insight. Even in the nonadiabatic model, any 
calculation of expectation values should be carried out 
within the adiabatic approximation. 



II. NONADIABATIC HEISENBERG MODEL 

A. General 

Consider a set of fi energy bands in a metal with the 
paramagnetic space group G, the paramagnetic group 



M p = G + KG 



(2.1) 



(with K denoting the operator of time inversion) , and \x 
atoms at the positions 



T = R + Qi 



(2.2) 



per unit cell, where R and Qi (i = 1 to fx) denote the 
vectors of the Bravais lattice and the positions of the zth 
atom within the unit cell, respectively. The energy bands 
of this set are assumed to belong to the narrowest partly 
filled bands of this metal (while it is not demanded that 
all the narrow, partly filled bands belong to it). 

Assume that the symmetry of the Bloch functions of 
the considered set of energy bands allows the construc- 
tion of either Wannier functions 



Wf m (r,t) = Wi(f '- R - Pi)u m {t) 



(2.3) 



symmetry-adapted to a magnetic subgroup M of M or 
spin dependent Wannier functions 



Wf m (r, t) = Wi m (r- R- pi, t) 



(2.4) 



symmetry- adapted to M p . The former functions arc 
defined appendix A, case (a), their symmetry is given 
in Eqs. (A17) and (A19), the latter are defined in ap- 
pendix A, case (b), and their symmetry is given in 
Eqs. (A25) and (A26). The functions u m (t) are Pauli's 
spin functions, see Eq. (4.7), t is the spin coordinate, and 
the (crystal) spin label m — ±i distinguished between 
the two functions at the same position T . In cither case, 
the Wannier functions form a unitary transformation of 
the exact Bloch functions of the considered set of /i en- 
ergy bands, are situated on the atoms (with the positions 
T), and are as well localized as possible. 



The three postulates of the nonadiabatic 
Heisenberg model 



Let be 



H = H 



HF 



H, 



ci, 



(2.5) 



the electronic Hamiltonian in the considered set of energy 
bands with Hhf and 



H cb = ^2(f 1 ,m 1 ;f 2 ,m 2 \H cb \f[,m / 1 ;f 2 ,m 2 ] 

f,m 



fimi f 2 m 2 ^2 m 2 7\' m i 



(2.6) 



representing the Hartree-Fock and Coulomb energy, re- 
spectively. The fermion operators c- and Cf m create 
and annihilate electrons with (crystal) spin m in the lo- 
calized states \T,m) represented by the Wannier func- 
tions Wf m (f, t). Other contributions to H from the elec- 
trons not belonging to the considered set of bands, are 
neglected even as are spin-orbit effects. 
Hcb may be written as 



Hcb = H c + H e 



H z 



(2.7) 



with the operator of Coulomb repulsion H c containing all 
the matrix elements of Hcb with 



fi = f[ and f 2 =f' 2 , 



(2.8) 



the exchange operator H ex containing the matrix ele- 
ments with 



7\ = f' 2 and f 2 = f[, 



(2.9) 



and H z comprising the remaining matrix elements, i.e., 
the matrix elements with 



{f u f 2 }^{fi,f^}. 



(2.10) 



The interaction H z is of great importance within the 
NHM. In order to discuss the effect of H z , consider the 
operator 



H' = H 



HF 



H r + H P 



(2.11) 



obtained from the complete Hamiltonian H by putting 
H z = 0. 

As is well-known, the Coulomb repulsion of two elec- 
trons occupying localized states at the same atom is 
larger than the Coulomb repulsion of two electrons at 
different atoms. 14 The electronic motion in the ground 
state \G') of H' has an "atomiclikc" character when the 
Coulomb repulsion between the localized states deter- 
mines the electronic motion in \G') to such an extend 
that the probability to find two electrons (with different 
spin directions) on the same atom is markedly smaller 
than in case of a purely bandlike motion. In this context, 
I speak of a "purely bandlike" motion when the proba- 
bility to find an electron in the localized state \Tm) is 
independent of whether or not the other state \T, —m) 
is occupied. In this case, the ground state consists of 
configurations with nearly random occupation. 

The atomic- or bandlikc character depends on the 
mean time of stay r ~ h/A of the electrons at the atoms 
and, hence, on the bandwidth A. For A — > we have 
t — > oo; the metal becomes a Mott insulator represent- 
ing a perfectly atomiclike state. For A — > oo, on the 
other hand, we have r — ► and, hence, a purely band- 
like character of the electrons. Thus, as is well-known, 
the electrons in partly filled energy bands tend to a more 
atomiclike behavior with decreasing bandwidth and to a 
more bandlike behavior with increasing bandwidth. 15 

Now assume the considered energy bands to be suffi- 
ciently narrow that the ground state \G') of H' clearly 
has atomiclikc character. The matrix elements of H z sat- 
isfy neither Eq. (2.8) nor (2.9). Thus, the interaction H z 
annihilates two electrons in localized states at the posi- 
tions T[ and T 2 and creates at least one of them at the 
new positions T\ or T 2 . Hence, unlike H c or H exi the op- 
erator H z generates transitions between adjacent local- 
ized states which lead to configurations with a more ran- 
dom occupation. Consequently, the probability to find 
two electrons at the same position T will be larger in the 
ground state \G) of the complete Hamiltonian 

H = H' + H Z 

[given in Eq. (2.5)] than in the ground state \G') of H' . 
Therefore, the total Coulomb repulsion energy in \G) is 



larger than in \G') and we may assume that in sufficiently 
narrow bands the ground state energy E of H is greater 
than the ground state energy E' of H' . 

Still "sufficiently" narrow means that the Coulomb re- 
pulsion between the localized states determines the elec- 
tronic motion in \G'). It is difficult to decide how narrow 
such sufficiently narrow bands should be. However, we 
know that there is strong theoretical and experimental 
evidence that, e.g., the d electrons of the transition met- 
als exhibit behavior of both the band and the Heisenberg 
model. 15 Therefore, we may suppose that E > E' is valid 
in the narrowest bands of the metals. This supposition 
leads to the first postulate of the NHM. 

Postulate 1 In the narrowest, partly filled energy bands 
of the metals the transitions generated by H z are ener- 
getically unfavorable, i.e., we have 



(G\H\G) > (G'\H'\G'), 



(2.12) 



where \G) and \G') denote the exact ground states of H 
and H' , respectively. 

The particular form of the matrix elements of H z shows 
that it represents a short-ranged interaction which cru- 
cially depends on the exact form of the localized func- 
tions. This fact suggests that only small changes of the 
localized electronic orbitals are required to prevent the 
transitions generated by H z . However, such modified or- 
bitals do not exist within the adiabatic approximation 
because such modifications yield localized charge distri- 
butions which are not symmetric at any moment. As 
a consequence, the nuclei become accelerated in varying 
directions. Hence, we replace the (adiabatic) localized 
states (represented by the Wannicr functions) by more 
realistic nonadiabatic localized states 



\T, m, v) 



(2.13) 



which take into account the motion of the nuclei. The 
new quantum number v labels different states of motion 
of the center of mass of the nucleus and the electron 
occupying the state \T, to, v). 12 

The nonadiabatic Hamiltonian H n may be written as 



rjn Tj | Tjn 

n — Hhf + ti Ch 



(2.14) 



where the Coulomb interaction now has the form 

H Cb = ^2(f 1 ,m 1 ,n;f 2 ,m2,n\Hcb\T[,m' 1 ,n;f2,m' 2 ,n) 

(2.15) 



T,m 

xc nf c" f c™ c" 

fimi f 2 m 2 T^m' 2 T[m' x ' 



The new fermion operators c™ and c™ create and an- 

Tra Tm 

nihilate electrons with crystal spin to [see appendix A, 
case (b)] in the nonadiabatic localized states \T, m, n). 



The matrix elements of H£ b are integrals 

(Ti, mi, n; T 2 , to 2 , n\Hcb\T[, m' 1 , n; T 2 , m 2 , n) 

= — ^2 / (fi,m 1 ,n\f,t,q)(f 2 ,m 2 ,n\r' ' ,t' ,q" ) 
tv •* 

x jzr Kn(r, t, g \T[, mi, n) (?', t', q'\T^m' 2 , n) 
\r — r | 

xdrdr'dqdq" (2.16) 

over nonadiabatic localized functions of the form 



(r,t,q\T,m,n), 



(2.17) 



where v — n labels the nonadiabatic states which satisfy 
the following Eq. (2.18), and the new coordinate q stands 
for that part of the motion of the center of mass of the 
localized state \T, m, n) which nonadiabatically follows 
the motion of the electron occupying this state. We may 
imagine that q denotes the acceleration of the nucleus 
(together with the core electrons). 

Within the nonadiabatic localized states \T, m, v) the 
electrons possess considerably more room to move than 
within the adiabatic states. While in the adiabatic ap- 
proximation the symmetry operators P(a) act on r and 
t alone, in the nonadiabatic system these operators act 
on r, t, and the acceleration q of the nuclei, see Eq. (B9). 
The nonadiabatic localized functions have no definite 
transformation properties under space group operations 
acting only on r and t. Hence, the symmetry of the adia- 
batic and nonadiabatic localized states [given by Eq. (Bl) 
for the fermion operators] may be interpreted as follows. 

Within the adiabatic system the electrons move on or- 
bitals being symmetric with respect to the lattice at any 
moment. Within the nonadiabatic system, on the other 
hand, the orbitals are still symmetric on the average of 
time, but not at any moment. This statement is indepen- 
dent of the absolute value \q\ of the acceleration of the 
nuclei, i.e., it is independent of whether or not the mass 
of the electrons is markedly smaller than the mass of the 
nuclei. 

Thus, the introduction of the new quantum number 
v allows the electrons to move in a potential depending 
on which of the adjacent localized states are occupied 
and on the present positions of these electrons. Hence, 
within the nonadiabatic system the electrons should be 
able to avoid the transitions generated by H z by an ap- 
propriately modified motion, if these transitions are en- 
ergetically unfavorable, i.e., if the relation (2.12) is true. 
Thus, as a consequence of relation (2.12), all the matrix 
elements of H^ b which neither satisfy Eq. (2.8) nor Eq. 
(2.9) should vanish. 

For this reason, we suppose that the transitions gener- 
ated by H z are artifacts of the adiabatic approximation 
and do not happen in the (true) nonadiabatic system if 
relation (2.12) is satisfied. This supposition leads to the 
second postulate of the NHM. 



Postulate 2 If relation (2.12) is true, the Coulomb in- 
teraction H^ b does not generate transitions between ad- 
jacent localized states, i.e., 

(Ti,mi,n;T2,m2,n\Hcb\T{,m' 1 ,n;T2,m' 2 ,n) = 

(2.18) 
for 

{fi,f 2 }^{f/,f^} 

and for special nonadiabatic localized functions 



(r, t, q \T, m, n) 



labeled by v = n. 



At the transition from the adiabatic to the nonadi- 
abatic system, the total energy of the electron system 
decreases by 



AE=(G\H\G)-(G'\H'\G') 



(2.19) 



if we neglect the energy of the nonadiabatic motion of the 
nuclei and the energy change caused by the slight modi- 
fication of the electronic orbitals within the nonadiabatic 
states. 

As a consequence of Eq. (2.18), the commutation 
properties of the operator H^, b depend on the symme- 
try of the nonadiabatic localized states. Since only small 
modifications of the adiabatic electronic orbitals are re- 
quired to prevent the transitions generated by H z , we 
can assume that the nonadiabatic Hamiltonian H n has 
the same commutation properties as the adiabatic Hamil- 
tonian H' given in Eq. (2.11). This is the third (and last) 
postulate of the NHM. 

Postulate 3 If relation (2.12) is true, the nonadiabatic 
Hamiltonian H n has the same commutation properties as 
the adiabatic Hamiltonian H' , i.e., 



[H',P] 



¥= 



o 



[H n ,P] 



+ 



0. 



(2.20) 



where P stands for any symmetry operator. 



As a consequence, the nonadiabatic localized functions 
have the same symmetry as the (adiabatic) Wannier func- 
tions Wrp (f, t) . 



III. SYMMETRY OF THE OPERATOR H' 

According to its definition, the operator H' arises from 
the complete adiabatic Hamiltonian H in Eq. (2.5) by 
putting 



H z = 0. 



(3.1) 



This equation does not state that H z is neglected, but 
that H z is put equal to zero. By this step, the commuta- 
tion properties of the operator H' depend on the symme- 
try of the Wannier functions, whereas the commutation 



properties of the complete adiabatic Hamiltonian H are 
independent of the symmetry of the used basis functions. 
The nonadiabatic matrix elements of H z , however, vanish 
within the NHM, see Eq. (2.18). 



Case (a): The Wannier functions are 
symmetry-adapted to a magnetic group 

If the Wannier functions are symmetry-adapted only 
to a magnetic subgroup M of the paramagnetic group 
M p , the symmetry of the operator H' is given by 



[H',P(a)]=0 for a e M 



(3.2) 



and 



[H', P(a)} ^ for a £ (M p - M), (3.3) 



where 

M p = G + KG 

stands for the paramagnetic group (and G is the space 
group). The symmetry operators P(a) are given in 
Eq. (A3) and K denotes the operator of time inversion. 
Especially, in this case (a) H' does not commute with K, 



[H',K}^0, 



(3.4) 



since K £ (M p - M). 

The first equation (3.2) is valid since the complete 
Hamiltonian H commutes with P(a) and also the oper- 
ator P(a)H'p- 1 (a) complies with Eqs. (2.8) and (2.9) if 
a £ M since the fermion operators P(a)c~ P^ 1 (a) and 
P(a)cf m P~ 1 (a) are for all the a £ M linear combina- 
tions of operators c- ; , and Cf,, ,, respectively, labeled 

by the same position T", see the equations (Bl) and (B2). 
Within the NHM it is important that H' does not com- 
mute with P{a) for a £ (M p — M). In the case (a) con- 
sidered in this section, the Wannier functions cannot be 
chosen in such a way that they are symmetry-adapted to 
a group M containing the operation a as well as all the 
elements of M. Consequently, for a £ (M p — M), the 
fermion operators P(a)c~ P^ 1 (a) do not comply with 
Eq. (Bl) or (B2), but are linear combinations 



P(a)cl P~ 



\a) 



E 

T'm' 
J 



1 ' TCI' .1 m^ t rpt 



T'm' 



(3.5) 



of at least two operators c- , with different labels T[ and 

Tg. We show that therefore the operator P(a)H'P^ 1 (a) 
has matrix elements violating Eq. (2.8) or (2.9). 

Consider a fermion operator combination belonging to 
the Coulomb interaction of H' , say 



O 



C T 1 C f 2 C '?2 C f 1 - 



(3.6) 



and assume for a special a £ (M p — M) the sum in 
Eq. (3.5) to consist of two summands, 



Pclp- 1 =o-ct + 6-ct, 



(3.7) 



labeled by the different positions U and V. In Eq. (3.7) 
we use the abbreviation P = P(a) and drop the index m 
since it does not matter here. With Eq. (3.7) we obtain 

pop- 1 = pd p- l PcL p- l p Cr? p- l p CrP p- 1 

Ti T 2 ^ 2 Ti 

x( a * Cc72 +6* Cl72 )( a * C[7i +&* Cl7i ) (3.8) 

where a^O and 6^0. For instance, the operator prod- 
uct 

C U! C V2 C U2 C Ul 

belonging to POP^ 1 and, hence, to PH'P -1 , violates 
Eq. (2.8) since Vi ^ U?.. Consequently, H' does not 
commute with P as expressed by Eq. (3.3). In the same 
way, H' does not commute with P when there are more 
than two summands on the right hand side of Eq. (3.7). 



Case (b): The Wannier functions are spin dependent 
and symmetry-adapted to the paramagnetic group 

If we consider spin dependent Wannier functions 
symmetry- adapted to the paramagnetic group M p , then 
we have 



[H', P(a)] = ior ae M 1 



and, especially, 



\H' K] = 0. 



(3.9) 



(3.10) 



However, in this case (b), the operator H' has matrix 
elements with 



mi + TO2 7^ TO i 



(3.11) 



because the coefficients f sm (q,k) in Eq. (A21) cannot 
be chosen independent of k, see appendix A, case (b). 
Therefore, H' does not conserve the crystal spin and, 
hence 16 , does not commute with the operators M(a) of 
the crystal spin defined in Eq. (A30), 



[H',M(a)]*0 

for at least one a £ Gm- 17 



(3.12) 



IV. SYMMETRY OF THE NONADIABATIC 
HAMILTONIAN H n 

A. Magnetic and paramagnetic group 

The nonadiabatic Hamiltonian H n has the same com- 
mutation properties as the adiabatic operator H' , see 



Eq. (2.20). However, the symmetry operators P(a) now 
act on r, t, and on the new coordinate q of the nonadia- 
batic localized functions, see Eq. (B9). 
Hence, we have 



(4.1) 
(4.2) 

(4.3) 

(4.4) 

(4.5) 



[H n ,P(a)}=0 foraeM, 

[H n ,P(a)} ^ for a £ (M p - M), 
and, especially, 

in the case (a) of the preceding section III, and 

[H n ,P(a)]=Q for a EM P , 
especially, 

[H n ,K} = 
in the case (b) of the preceding section III. 

B. Crystal spin 



The nonadiabatic fermion operators in Eq. (2.14) are 
no longer labeled by the spin quantum number s. Hence, 
within the nonadiabatic system, the exact Fermi excita- 
tions are no longer purely electronic states but localized 
states of well-defined symmetry which are occupied by 
electrons carrying with them some nonadiabatic motion 
of the nuclei. 

Let be S(a) with 

S(a)u s (t) = u s (o~ t) = \ d s i s (a)u s i(t) for a £ 0(3) 

(4.6) 
the operators turning the electron spin, where the func- 
tions 



u s (t) 



(4.7) 



are Pauli's spin functions with the spin quantum number 
s = ±2 and the spin coordinate t = ±5, and the matrices 
[d s 'g(o;)] are the representatives of the two-dimensional 
double-valued representation D 1 i 2 of the three-dimensi- 
onal rotation group 0(3). 

The adiabatic Hamiltonian H given in Eq. (2.5) com- 
mutes with the operators S(a), 



[H,S(a)]=0 for a e 0(3). 



(4.8) 



This equation expresses the conservation law of the spin 
angular momentum within the adiabatic system. 

The nonadiabatic Hamiltonian H n , on the other hand, 
does not commute with the operators S(a) (for a ^ E) 
since the nonadiabatic fermion operators arc no longer 
labeled by the spin quantum number s. Hence, as a con- 
sequence of the (small) shift of the Fermi character at the 



transition from the adiabatic to the nonadiabatic system, 
the electron spin angular momentum is no longer a con- 
served quantity. Now there exists an interaction between 
the electron spins and the nonadiabatic motion of the 
nuclei. 

However, even in the nonadiabatic system there should 
exist a conserved quantity related to the conservation law 
of angular momentum. Thus, the equation (4.8) should 
be replaced by an analogous equation 



[H n ,M(a)}=0 for a e G 



M 



(4.9) 



in the nonadiabatic system, where TL n stands for the com- 
plete nonadiabatic Hamiltonian. The group Gm and the 
operators M(a) are defined in Eqs. (A29) and (A30). 
They act on the quantum number m of the nonadiabatic 
localized states \T, to, n) in the same manner as the oper- 
ators S(a) act on the spin quantum number s of Pauli's 
spin functions u s (t), cf. Eq. (B3). Therefore, these oper- 
ators may be called the symmetry operators of the "crys- 
tal spin" and to may be called the quantum number of 
the crystal spin. This is in analogy to the wave vector k 
of the Bloch functions which is sometimes referred to as 
"crystal momentum" in order to distinguish it from the 
momentum p. 

In the case (b), i.e., if we consider spin dependent Wan- 
nier functions symmetry-adapted to the paramagnetic 
group M p , the adiabatic operator H' does not commute 
with all the operators M(a), see Eq. (3.12). Hence, also 
the nonadiabatic Hamiltonian H n as defined in Eq. (2.14) 
does not conserve the crystal spin, 



[H n ,M(a)}^0 



(4.10) 



for at least one a E Gm- 17 It is one of the most interest- 
ing problems of the NHM to interpret this equation, see 
Sec. VB, case (b). 



DISCUSSION 



A. Crystal electrons 



The nonadiabatic Hamiltonian H n , on the other hand, 
has an important feature which distinguishes it from any 
adiabatic Hamiltonian H: the commutation properties 
and the spin dependence of H n depend on the symmetry 
and spin dependence of the nonadiabatic localized func- 
tions. This is because the nonadiabatic localized func- 
tions have a physical meaning going beyond the meaning 
of pure basis functions: they represent states that are 
really occupied by the electrons in the way described by 
Mott 14 and Hubbard 15 : the electrons occupy the nona- 
diabatic localized states as long as possible and perform 
their band motion by hopping from one atom to another. 
Such a band motion is generally referred to as atomiclikc 
motion. 

Within the NHM we may extend this picture of the 
atomiclikc electron. Here the whole localized state 
\T, to, n) behaves like a moving particle, say "crystal elec- 
tron" , with the local coordinate T and the crystal spin to. 
The spatial extend of the crystal electron is determined 
by the charge distribution of the localized state and the 
crystal spin is a conserved quantity. 

First consider the picture of the crystal electron within 
the adiabatic approximation. Both operators H c and H ex 
[given in Eq. (2.7)] represent interactions between crystal 
electrons and, hence, are in accordance with this pic- 
ture. The interaction H Zl on the other hand, contradicts 
the picture of a moving crystal electron because it de- 
stroys these new particles. Since H z is a short-ranged 
interaction, we may say that within the adiabatic system 
"the crystal electrons become destroyed at the slightest 
touch" . 

Within the NHM, on the other hand, Eq. (2.18) is 
valid. The crystal electrons are stable in the nonadiabatic 
system because the Coulomb interaction does not gen- 
erate transitions between adjacent localized states. We 
may interpret Eq. (2.18) by stating that "the crystal elec- 
trons become slightly deformed but not destroyed at a 
touch" . The crystal electrons now have a certain elastic- 
ity protecting them from being destroyed at any collision. 
In this context, the stability of crystal electrons increases 
with decreasing band width. 



The NHM has been developed in order to interpret the 
symmetry and spin dependence of the Wannier functions 
in metals. These Wannier functions form an exactly uni- 
tary transformation of the Bloch functions of a set of 
partly filled energy bands in the band structure of the 
metal of interest. 

In the framework of the adiabatic approximation, 
Wannier functions form nothing but a unitary basis 
within the considered bands. Hence, the commutation 
properties of the adiabatic Hamiltonian H are indepen- 
dent of the symmetry of the Wannier functions, since the 
symmetry of any adiabatic Hamiltonian is independent 
of the symmetry of the used basis functions. The symme- 
try and localization of the Wannier functions simplifies 
the calculation of the matrix elements of H , but has no 
further physical meaning. 



B. Outlook 

The symmetry of the nonadiabatic Hamiltonian H n 
is given in Sec. IV for two interesting cases (a) and (b) 
which shall be considered separately. 



Case (a): The Wannier functions are 
symmetry-adapted to a magnetic group 

Sets of narrow, roughly half-filled energy bands with 
Wannier functions symmetry-adapted to a magnetic 
group M as defined in the appendix A, case (a), have 
already been identified in the paramagnetic band struc- 
tures of iron 9 and chromium 8 . In both metals, this set 



consists of one "magnetic" band. In iron the related mag- 
netic group 



M = IA/mm'm' = C\ h + K{C 2x \Q}Cl h 



(5.1) 



is the group of the ferromagnetic state and in chromium 



M = P/4/mnc = D\ h + K{E\t}D, 



■a. 



(5.2) 



is the group of the commensurate spin-density-wave 
state. 

In the case (a) considered in this section, the symme- 
try of the nonadiabatic Hamiltonian H n is given by the 
Eqs. (4.1), (4.2), and (4.3). Equation (4.3) shows that 
H n does not commute with the operator K of time in- 
version. Therefore, H n cannot have a paramagnetic or 
superconducting ground state, since both states are in- 
variant with respect to the time inversion. 

Thus, the electrons of the magnetic band may gain the 
energy AE [given in Eq. (2.19)] only if the electron spins 
form a structure with the magnetic group M. This fact 
may be interpreted as follows. 8 ' 9 

The electrons of the magnetic band activate a spin de- 
pendent exchange mechanism producing a spin structure 
with the space group M. This is possible since, first, the 
electrons can modify their orbitals in the nonadiabatic 
localized states and, secondly, exchange integrals depend 
very sensitively on the exact form of the electronic or- 
bitals. Hence, the electrons of the magnetic band modify 
their orbitals in such a way that the exchange energy E ex 
is maximum for a spin structure with the group M. 

The condensation energy Ef, i.e., the energy difference 
between the paramagnetic and the magnetic state, is no 
longer given by the exchange energy E ex alone, but by 



Ef — AE + E ex . 
Hence, Ef may be positive even if E ex is negative. 



(5.3) 



Case (b): The Wannier functions are spin dependent 
and symmetry-adapted to the paramagnetic group 

Sets of narrow, roughly half-filled energy bands with 
spin dependent Wannier functions symmetry-adapted to 
the paramagnetic group M p as defined in the appendix 
A, case (b), have already been identified in the band 
structures of a great number of superconductors. 7,10 ' 11 
It is remarkable that such "superconducting" bands 
cannot be found in those metals (such as Li, Na, K, 
Rb, Cs, Ca Cu, Ag, and Au) which do not become 
superconducting. 11 

In the case (b) considered in this section, the symme- 
try of the nonadiabatic Hamiltonian H n is given by the 
Eqs. (4.4), (4.5), and (4.10). Eqs. (4.4) and (4.5) show 
that the ground state of H n has the correct symmetry of 
the paramagnetic group M p . Especially, H n commutes 
with the operator K of time inversion. Therefore, H n 
may have a paramagnetic or a superconducting ground 
state, but has not a magnetic ground state. 



Equation (4.10) shows that H n does not conserve the 
crystal spin angular momentum. Hence, the electrons of 
the considered bands may gain the energy AE [given in 
Eq. (2.19)] only if they couple to other excitations in such 
a way that the conservation of the crystal spin angular 
momentum is satisfied in the nonadiabatic system. This 
fact may be interpreted as follows. 7 ' 18 ' 19 

In isotropic materials, the electrons of the considered 
bands couple to the phonons. This is possible since, first, 
the symmetry of localized acoustic phonons shows that 
they are able to carry crystal spin angular momentum, 
secondly, the electron spins are coupled to the phonons 
via the nonadiabatic motion of the nuclei, and, thirdly, 
the resulting nonadiabatic Hamiltonian complies with the 
conservation law of crystal spin angular momentum. 

In anisotropic materials (consisting of one- or two- 
dimensional sublattices), phonons are not able to trans- 
port crystal spin angular momenta through the crystal. 
Here the electrons of the considered bands are forced to 
couple to energetically higher-lying boson excitations. 

At zero temperature, this spin-boson interaction con- 
strains the electrons of the considered bands in a new 
way to form Cooper pairs because the conservation of 
spin angular momentum would be violated in any normal 
conducting state. Apart from this participation of the 
conservation of spin angular momentum, the mechanism 
of Cooper pair formation within the NHM is identical to 
the familiar mechanism presented within the Bardeen- 
Cooper-Schrieffer (BCS) theory. 20 

The participation of the conservation law of spin an- 
gular momentum may be interpreted in terms of quan- 
tum mechanical constraining forces which constrain the 
electrons to form Cooper pairs. There is evidence that 
these constraining forces are necessary for the Hamil- 
tonian to have eigenstates in which the electrons form 
Cooper pairs. 

If this is true, then the BCS theory of superconductiv- 
ity is only applicable to superconducting bands as defined 
in the appendix A, case (b). When it is applied to other 
bands, it does not yield the absolute energy minimum in 
the Hilbert space. 
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APPENDIX A: SYMMETRY-ADAPTED 
WANNIER FUNCTIONS 

Consider a metal with the space group G and the point 
jroup Gq. The elements 



{a\t} 



(Al) 



of G consist of a point group operation a and a transla- 
tion vector 



t = f(a) + R 



(A2) 



which is the sum of the nonprimitivc translation r(a) as- 
sociated with a and a translation vector R of the Bravais 
lattice. 

The operators P(a) act on a wave function /(r, t) de- 
pending on the position r and the spin coordinate t ac- 
cording to 



P(a)f(f, t) = /(a _1 f - oT x t, aT x t), 

where the symbol aT x t is defined in Eq. (4.6). 
The effect of K is given by the equations 21 



(A3) 



Kf(f) - f*(r), (A4) 

where f(r) stands for any function of position, and 

Ku s (t) = g s u- s (t), (A5) 

with 22 

#±1/2 = T«- (A6) 

Case (a): Wannier functions symmetry-adapted to a 
magnetic group 

Consider a set of p energy bands in the paramagnetic 
band structure of a metal with p atoms per unit cell. The 
positions of the atoms are still written as 

f = R + g t , 

where R and Qi (i = 1 to p) denote the vectors of the 
Bravais lattice and the positions of the zth atom within 
the unit cell, respectively. 

Further, consider the magnetic group 

M = H + K{j\f(j)}H (A7) 

where H is a subgroup of G, 

H cG, 

K denotes the operator of time inversion, and {7|t(7)} 
is a space group element of G — H . 

Assume degeneracies to exist between the bands of the 
considered set of p energy bands and the bands not be- 
longing to this set. This assumption is always true since, 
in any metal, there arc degeneracies between the bands 
of any selected set of energy bands and the bands outside 
of this set. These degeneracies are caused by symmetry 
and may occur at points and lines of symmetry of the 
Brillouin zone. Further, assume 

• these degeneracies to be removed in the subgroup 
H of G (i.e., when the representations of G are 
replaced by the subduced representations of H) ; 



• the symmetry operation i^{7|r(7)} not to produce 
extra degeneracies between the bands of the con- 
sidered set and bands outside of this set; 

• unitary matrices S(K) [as defined in Eq. (4.16) of 
Ref. 2] to exist which satisfy the equations (4.17) 
and (4.28) of Ref. 2 and Eq.' (7.1) of Ref. 13; and 

• the positions pi of the Wannier functions [which are 
determined by these matrices S(K)] to be identical 
with the positions of the atoms. 

Then the coefficients gi q (k) in Eq. (1.1) may be chosen 
in such a way that the Wannier functions comply with 
the four conditions following Eq. (1.1) with the exception 
that they are no longer symmetry-adapted to the param- 
agnetic space group G but only to the subgroup H of G. 
That means, Eq. (1.3) is satisfied only for the elements 
a of the point group Hq of H. In addition to Eq. (1.3) 
we have 13 

Kw l ( 1 - l {r-R-p^ i )) =J2D ji (K^r) Wj (f-R-p i ) (A8) 
i=i 

where the matrix [Dji(K"f)] is the representative of Kj 
in the corepresentation of the point group 



M = H + K-fH 



(A9) 



of M which is derived from the representation Do of Hq 
in Eq. (1.3). 

Since there is exactly one Wannier function at each 
atom, the Wannier functions may be labeled by the po- 
sitions T of the atoms, 



Wf{r) = u>i(r- R - p t ), 



(A10) 



and the equations (1.3) and (A8) may be considerably 
simplified. 

Applying on both sides of Eq. (1.3) the operation a on 
r — R — pi we obtain 



a; 



{r-R-pJ^^DjiiaJwjiair-R-Pi)), (All) 



and the application of the operator P(a) [given in 
Eq. (A3)] on both sides of this equation (All) yields the 
equation 

P(a)wi(r — R — pi) = y Dji(a)wj(f — t — aR — api) 

i=i 

(A12) 
which applies to all the elements a £ H. 

As shown in Ref. 6, Eq. (A12) may also be written in 
the form 

P(a)iVi(r — R— pi) = y^ Dji(a)iVj (r — aR — pj — Rj(a)) 

i=i 

(A13) 
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with Rj (a) being translations of the Bravais lattice (de- 
pending on j and a), see Eq. (2.13) of Ref. 6. 

Comparing the Eqs. (A12) and (A13), we see that on 
the right hand side of Eq. (A13) there are only Wannier 
functions related to an atomic position pj within the unit 
cell for which the translation 



Rj (a) = api + t — Pj 



(A14) 



is a translation vector of the Bravais lattice. This can- 
not be true for more than one vector pj since all the pj 
are different and lie within the unit cell. It is true for ex- 
actly one pj because Rj (a) is a translation of the Bravais 
lattice if we put 



Pj = api + t 



(A15) 



and R + api + 1 is the position of an atom since it can be 
generated by the application of the space group operation 
{a\t\ on the atomic position a^ 1 R + pj. 

Consequently, on the right hand side of Eq. (A12) there 
is only one Wannier function, namely that function re- 
lated to the atom at position api + 1 within the unit cell. 
Hence, the sum on the right hand side of Eq. (A12) and 
(analogously) of Eq. (A8) consists of one summand only. 
The matrices [Dji(a)] in these equations have only one 
non- vanishing element, say dji(a), in each column which 
satisfies the equation 



\dji(a)\ = 1 

since the matrix [Dji(a)\ is unitary. 
Hence, Eq. (A12) may be written as 



P(a)Wf(r) — df(a)wf,(r) 



for a € H 



(A16) 



(A17) 



with 



T' = aT + t (A18) 

and Eq. (A8) yields 

KP{g)w f {r) = df{Ki)w fl {r) (A19) 

with 



and 



g = {7^(7)} 



f^^f+Tij), 



where the coefficients df(a) and df(Kj) have the abso- 
lute value 1, 



\df(a)\ = \d f (K 7 )\ = 1. 



(A20) 



It should be noted that the time inversion K does not 
belong to M. Therefore, it is not possible to choose the 
coefficients gi q {k) in Eq. (1.1) in such a way that the Wan- 
nier functions satisfy an equation analogous to Eq. (A19) 
by application of the time inversion operator K alone. 



Case (b): Spin dependent Wannier functions 
symmetry-adapted to the paramagnetic group 

a. Symmetry operators 

If in Eq. (1.1) we replace the Bloch functions ip? (r) 
by Bloch functions 



^4 



kqm 



(f,t)= Y, fsm(q,k)u s (t)ip Sq {f) (A21) 



with k dependent spin directions, we get "spin dependent 
Wannier functions" 



w im (r-R-pi,t) 
1 






V h 9=1 



■ <E( * + *Wfc)^, m (?,t), 



k 1= 



(A22) 



which are labeled by the additional quantum number 
m = ±5 of the crystal spin. The functions u s (t) de- 
note Pauli's spin functions as given in Eq. (4.7) and the 
coefficients f sm {q,k) form (for each k and q) a unitary 
two-dimensional matrix f(q, fc), 



If we have 



r\q,k) = f(q,k). 



fsm(q,k) = 5 S 



(A23) 



(A24) 



the two functions 4>t (r,t) (with m = ±|) are usual 
Bloch functions with the spins lying in +z and — z direc- 
tion, respectively. Otherwise, the functions <p^ m (r, t) 
still are usual Bloch functions with antiparallcl spins 
which, however, no longer lie in ±z direction. 

As in the preceding case (a) , consider a set of p energy 
bands in the paramagnetic band structure of a metal with 
p atoms per unit cell. 

The paramagnetic group M p of the metal may be writ- 
ten as 



M l 



G + KG 



where K still denotes the operator of time inversion. As- 



• the symmetry degeneracies between the bands be- 
longing the considered set and bands not belong- 
ing to this set to be removed when the single- 
valued representations of G are replaced by the cor- 
responding double-valued representations; 

• the time inversion symmetry not to produce extra 
degeneracies between the bands of the considered 
set and bands outside of this set; 
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• unitary matrices S(K) to exist which satisfy the 
equations (4.16), (4.17), and (4.28) of Rcf. 2 and 
Eq. (7.1) of Ref. 13, when the single-valued repre- 
sentations in these equations are replaced by the 
corresponding double-valued representations; and 

• the positions pi of the Wannicr functions to be iden- 
tical with the positions of the atoms. 

Then the coefficients f sm (q t k) and gi q (k) in Eqs. (A21) 
and (A22) may be chosen in such a way that also the 
spin dependent Wannier functions comply with the four 
conditions following Eq. (1.1). However, the coefficients 
fsm(q,k) cannot be chosen independent of k since the 
considered set of energy bands is isolated not before the 
single-valued representations of the Bloch functions are 
replaced by the related double- valued representations. [If 
the f sm (q,k) are independent of k, the Wannier func- 
tions in Eq. (A22) are usual Wannier functions which 
may comply with the four conditions following Eq. (1.1) 
only if the considered set of energy bands is already iso- 
lated when the Bloch functions are labeled by the single- 
valued representations of G.\ As an important conse- 
quence, the operator H' does not conserve the crystal 
spin, see Eq. (3.12). 

The symmetry of the spin dependent Wannier func- 
tions may be derived from the equations in Rcfs. 2, 6 
and 13 in the same way as we have derived the symme- 
try of the magnetic Wannier functions in the preceding 
case (a). We now get the equations 

+h 
p { a ) w f m ^^) = df(a) ^2 dm'm(a)Wf, ml (f,t) 

(A25) 
for a G G, and 

+ k 
Kw fm {r,t)=d f (K) J2 d m/m (K)w fm ,(f,t), (A26) 

where 



Wfjf.t) 



,(f- R- Pi,t) 



and 



f' = aT + t. 



The operators P(a) now act on f and t, see Eq. (A3), 
the matrices [d m i m (a)] are the representatives of the 
two-dimensional double-valued representation ~D\ii of 
the three-dimensional rotation group 0(3), the matrix 
[d m ' m (K)] is given by 



[d m >m{K)\ 



1 
-1 



(A27) 



[see, e.g., Table 7.15 of Ref. 21], and the c-numbers df(a) 
and df{K) still have the absolute value 1, 



The equations (A25) and (A26) are already given in 
Rcf. 10. Further, in appendix B of Ref. 10 simple equa- 
tions are given to identify sets of energy bands complying 
with all the condition given above. 



b. Operators of the crystal spin 

Define the "group of the positions p;" Gm to consist 
of all the a € Gq which satisfy the equation 



ap % + f (a) = pi + Ri 



(A29) 



for each pi, where Ri denotes a translation vector of the 
Bravais lattice, and define for all a € Gm symmetry 
operators of the "crystal spin" 

M{a) = P{{E\R- Ri})P[{a\T{a)}]P({E\ - R}) (A30) 

which depend on the position 

f=R + p t 

of the (spin dependent) Wannicr function on which they 
are acting. 

From Eqs. (A30) and (A25) we obtain the equation 

H 
M(a)wf m (f,t) = df(a) ^ d m > m (a)wf m ,(r,t) 
m '— — \ 

(A31) 
for a E Gm, showing that the operators M(a) leave 
unchanged the positions of the spin dependent Wannicr 
functions. 



APPENDIX B: SYMMETRY OF THE FERMION 
OPERATORS 

1. Adiabatic fermion operators 

The adiabatic fermion operators c- and Cf m create 
and annihilate electrons in localized states represented 
by the the Wannier functions Wf{r)u m (t) [in the case (a) 
of magnetic Wannier functions] or Wf (r,t) [in the case 
(b) of spin dependent Wannier functions]. Hence, their 
symmetry is determined by the equations (4.6), (A17), 
(A19), (A25), and (A26). From these equations we get 



+ * 



Pia^p-^a) 



with 



\df(a)\ = \d f (K)\ = 1. 



(A28) 



df{a) Yl 



T' = aT + t 



.(«)«, (Bl) 
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and 



and in Eq. (B2) we have 



KP(g)c^[KP(g)]-'=d f (K 7 ) ]T d m , m (Kj)cf>l, 

(B2) 



g = {£|0} 



(B7) 



with 



?'=-,? +7^), 



where the superscript (n) of the fermion operators should 
be disregarded in this section. In addition, the operators 
of the crystal spin satisfy the equation 



M(a)c^M-\a) = d f (a) ]T <W«)cKJ (B3) 

m'=-i 

for all the elements a of the group of the positions Gm , 
sec Eqs. (A29) and (A30). 

The coefficients df(a) and df(K'y) still have the ab- 
solute value 1, 

I4r(«)l = \ d f (#7)1 = 1, 

and the matrices [d m > m (a)] and [d m > m (Kj)] still are rep- 
resentatives of the two-dimensional double- valued repre- 
sentation D1/2 of the three-dimensional rotation group 
0(3) and the corepresentation of 0(3) + KO(3) derived 
from -D1/2, respectively. 

In the case (a) of the magnetic Wannier functions, 
Eq. (Bl) is valid for all the elements a in H, 



ae H, 
and in Eq. (B2) we have 

g = (7|f(7)}. 



(B4) 



(B5) 



In the case (b) of the spin dependent Wannier functions 
Eq. (Bl) is valid for all the elements a in G, 



where E denotes the identity element of Go- In the latter 
case the matrix [d m ' m (K)] is given in Eq. (A27). 



2. Nonadiabatic fermion operators 

Within the NHM, the Wannier functions Wf(r)u m (t) 
[in the case (a) of magnetic Wannier functions] or 
Wf m (f,t) [in the case (b) of spin dependent Wannier 
functions] are replaced by nonadiabatic localized func- 
tions, 



(f,t,q\T,m,n), (B8) 



ae G, 



(B6) 



Wf(r)u m (t) 

Wfjr,t) 

having the same symmetry as the Wannier functions. 
However, in the case of the nonadiabatic localized states, 
the symmetry operators P(a) act on r, t, and on the new 
coordinate q according to 

P(a) (r, t, q \T, m, n) 

= (a - f — oT t,a~ t,a~ q\T,m,n), (B9) 

where the symbol a~ 1 t is defined in Eq. (4.6), and the 
application of K yields 

K(r, t, q \T, m, n) = g m (r,t,q\T,-m,n)*, (BIO) 

with g m being given in Eq. (A6). 

With these redefinitions of the symmetry operators, 
the symmetry of the nonadiabatic fermion operators 
c™ is also given by the equations of the preceding ap- 
pendix B 1. The superscript (n) of the fermion operators 
in these equations shall indicate that they are valid for 
both the adiabatic and nonadiabatic fermion operators. 
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